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In this mini-review, we show that a lot of theoretical efforts have been made for the theoretical 

study of two body hadronic _B( 3 j and B c decays. In addition to many next-to-leading order or even 

next-to-next-to leading order a s corrections made, we also study many of the previously unknown 

next-to-leading order power corrections. While the former corrections are theoretically solid, the 

latter corrections are phenomenologically more important. In the QCD factorization approach based 

on collinear factorization, there is difficulty to deal with the power correction diagrams due to the 
i — i 

endpoint singularity. Thus many of these analysis use phenomenological method. In the perturbative 
£S) , QCD approach based on kr factorization, the endpoint singularity is killed by including the quark 

transverse momentum. Therefore we can calculate the annihilation type diagrams quantitatively, 
which give the right sign for the direct CP asymmetry parameters. More and more _B( S ) decays 
1 channels, especially the pure annihilation type decays have been measured by the recent experiments 

^ ■ to confirm our theoretical predictions. More channels are predicted for future experiments, such as 

\ the charmless and charmed B s and B c hadronic decays and the decays involving a scalar, axial 

vector, even a tensor meson in the final states. 

I. INTRODUCTION 

In the LHC era of particle physics, heavy flavor physics is still one of the hot topics in particle physics. Most of 
the standard model free parameters belong to the flavor part. One of the four experiments, the LHCb experiment 
focus mainly on the B physics, although other two experiments-ATLAS and CMS also providing us rich data of B 
" physics. In fact, while the two B factories continue to analyze data, two new super B factories are preparing for 100 
i more luminosity. With the big achievements in experiments, the theoretical improvement of B physics has been a 
' little bit slow down. 



d 



Although many of the next-to-leading order or even next-to-next-to-leading order a s corrections have been done in 
the perturbative QCD approach [l|, the QCD factorization approach and the soft-collinear-effective theory study Q, 
the phenomenologically more important power corrections have not been systematically studied in these approaches. 
It has been demonstrated that in many non-leptonic B decays, the previously missed power correction, such as 
annihilationtype diagrams are very important in the direct CP asymmetry and polarization study of vector meson 
final states [3j. As pointed out, these power corrections are not calculable in the QCD factorization approach with 
endpoint singularity In fact, the well defined perturbative QCD approach with kx factorization Q can calculate 
all of these kinds of diagrams without ambiguity. 

In two body hadronic Bt s -\ decays, the light final state mesons and their constituent quarks inside are collinear 
objects at the rest frame of B/ s \ meson. The light spectator quark in the i?( s ) meson is rather soft. Therefore a 
hard gluon is needed to transform it into a collinear object. The dominant contribution is perturbative. However, 
the calculation is sometimes divergent at the endpoint of the meson distribution amplitudes. To deal with this 
unphysical singularity, one leaves the form factor diagrams to non-perturbative contributions in QCD factorization 
approach and soft-collinear effective theory (SCET)[6|. For the annihilation type diagrams, the QCD factorization 
approach just parameterizes it as free parameters; while SCET argues it as small power corrections. In fact, the 
quark carries very little longitudinal momentum at the endpoint, therefore the transverse momentum of quark is no 
longer negligible. In the perturbative QCD approach, we keep the transverse momentum of quark, which acts as a 
natural regulator of the endpoint divergence. Including another momentum scale (transverse momentum) in QCD 
will produce large double logarithms in the perturbative calculations, so that we have to use the renormalization 
group equation to do the resummation. A Sudakov factor is produced, which suppresses the endpoint contributions 
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to keep the perturbative calculation healthy. 

In this paper, we will show that many of the leading order and also next-to-leading order perturbative QCD 
calculations of Br 8 \ and B c decays have been tested in the experiments and more and more channels are predicted 
for future experiments. All of the charmless or charmed meson decays of B s meson have been studied some years 
ago in the perturbative QCD approach [7, 8] and the QCD factorization approach 0, Some of the charmless B s 
decays are also discussed in the soft-collinear effective theory 0, [nj ■ Recently, a scalar meson, axial vector meson, 
or tensor meson involved in the final states are also studied. Since the LHC experiment will produce numerous B c 
mesons, the intensive study of B c meson hadronic decays are also studied, which include the charmless final states, 
one or two D mesons in the final states and even final states involving scalar or tensor mesons. 

II. FACTORIZATION METHOD IN QCD FOR HADRONIC B (s) DECAYS 

All the hadronic B/ s \ decays are weak decays. In the quark level, the weak transitions are summarized as effective 



four quark operators with QCD corrections. The weak effective Hamiltonian Heff can be written as 11 1 

G F 



Heff = 



\ E KbVix [Ci(/i)o?( M ) + c 2 ( M )o 2 V)] - vsytx 



i=3 



(1) 



with V q b(x) an d V t b(x) (X = d, s) the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The current-current 
(tree) local four-quark operators are: 

Of = {b a qp)v-A{qf3X a ) v ^ A , 0% = (b a q a )v-A(qpXp)v-A] (2) 

the QCD penguin operators are 

3 = (baX^v^A^iq'pq'^v^A, 4 = (b a X l3 )v-A^2(cfpq' a )v-A, (3) 

q' q> 

5 = (baX^v^A^iq'pq'pjv+A, O e = (b a Xp) V -A^2%q' a )v+A, (4) 

and the electro-weak penguin operators are 

O7 = ^(b a X a ) v ^ A ^2eq l (q' q' l3 ) v+A , O s = -(b a Xp) V -A ^ e <?' {<?0<l'a)v+A, (5) 

q> q' 

09 = -Ma)v-A^e,'(^g|i)v-i, O10 = -{b a Xp) V - A '^2 i eq,(q'pq' a ) V -A, (6) 

q' q' 

where a and /? are the color indices and q' are the active quarks at the scale mt, i. e. q' — (u, d, s, c, b). The left-handed 
and right-handed currents are defined as (b a qp)v-A = &q7^(1 — lb)qfi and (q'rQ'oJv+A = 9^7^(1 + 75)?qi respectively. 
The Wilson coefficients C-s are calculated by the renormalization group equations to include the next-to-leading 
order QCD corrections. 

When doing the hadronic matrix elements calculations, we need to deal with scales around ^/^QCD™6, which 
is usually called factorization scale. Scales below this scale are categorized as non-perturbative physics, which 
is described by form factors or meson wave functions. In the generalized factorization approach [12j . the QCD 
corrected Wilson coefficients are scale dependent while the matrix elements described by the form factors are scale 
independent physical quantities. It is later improved by the QCD factorization approach Q to remove the scale 
dependence in the decay amplitudes by including the vertex corrections together with the hard scattering diagrams 
shown at FigOJc) and (d). In this collinear factorization based theory, however, there is endpoint singularity at the 
higher twist calculation and the annihilation type diagrams shown in Fi^2] Those annihilation type diagrams are 
later proved to be important Similar to the QCD factorization approach, the soft-collinear effective theory 0] 
also leave part of the soft contribution in the form factor diagrams shown in FigJIJa) and (b) as non-perturbative 
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FIG. 1: The leading order emission tip Feynman diagrams 
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FIG. 2: The leading order annihilation type Feynman diagrams 

contribution. This makes the soft-collinear effective theory less predictive, since it requires more free parameters to 
be determined by experiments 0, [l(3] • 

In the perturbative QCD approach, we keep the transverse momentum of quarks to kill the endpoint divergence. 
Therefore, we have one more scale, i.e. the quark transverse momentum than the QCD factorization approach and 
SCET. The factorization formula in pQCD approach is then 



C(t) <g> H(x, t) ® $(x) ® exp[-s(P, b) 



—lq(a s {p))] 

1/6 M 



(7) 



where C(t) are the corresponding Wilson coefhcients of four quark operators, which include the dynamics of physics 
larger than nib scale. The Sudakov evolution exp[— s(P, b)} [13J are from the resummation of double logarithms 
ln 2 (P6), with P denoting the dominant light-cone component of meson momentum. ^ q = —a s /n is the quark 
anomalous dimension in axial gauge. All non-perturbative components are organized in the form of hadron wave 
functions which may be extracted from experimental data or other non-perturbative method, such as QCD sum 

rules. Since non-perturbative dynamics has been factored out, one can evaluate all possible Feynman diagrams for the 
six-quark amplitude straightforwardly, which include both factorizable and non-factorizable emission contributions 
shown in FigfT] Factorizable and non-factorizable annihilation type diagrams shown in Figj2] are also calculable 
without endpoint singularity. 

A. Wave Functions of the B s Meson 



In order to calculate the analytic formulas of the decay amplitudes, we need the light cone wave functions decom- 
posed in terms of the spin structure. In general, the light cone wave functions are decomposed into 16 independent 
components, l a p, j£p, cr^, (7^75)0,3, j5a/3- If the considered meson is the heavy B or B s meson, the meson 
light-cone matrix element can be decomposed as [lil . Il5| 



d 4 ze lk ^(0\b a (0)s^z)\B s (P B J) 

u 



\/6 I 



4>B s {h) - ^ <t>B,(ki] 



(8) 



Pa 
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FIG. 3: B a meson distribution amplitudes. The solid-, dashed-, and tiny-dashed- lines correspond to ojb = 0.45 GeV, 
0.5 GeV, and 0.55 GeV, respectively. 



From the above equation, one can see that there are two Lorentz structures in the B s meson distribution amplitudes, 
that obey the following normalization conditions 



d 4 h 



fB. 



4>b s (fa 



0. 



(9) 



In general, one should consider these two Lorentz structures in the calculations of B s meson decays. Since the 
contribution of <j) Bs is numerically small |l6j |. its contribution can be neglected. With this approximation, we only 
keep the first term in the square bracket from the full Lorentz structure in Eq. (|8]) 



= —(P Bs +M B .)75& J ,(A;i). 
V6 



(10) 



Usually the hard part is always independent of one of the k^ and/or k x . The B s meson wave function is then a 
function of the variables k^ (or k^) and k^ only, 



4>B a (ki,K) = J ^L0 Ss (fc+,fcr,fc r ) 
The B s meson's wave function in the 6-space can be expressed by 



®b s (x,b) 



V6 



[Pb b 75 + M Bs 75 ] <f) Bs (x,b), 



(11) 



(12) 



where b is the conjugate space coordinate of the transverse momentum k . 

In this study, we use the model function similar to that of the B meson which is 

„2 



4> Bs (x,b) = N Bs x (1 - x) exp 



Ml 



1 



2ul 



{uj b by 



(13) 



with N Bs the normalization factor. In recent years, a lot of studies have been performed for the B^ and B^ decays 
in the pQCD approach [5]. The parameter uib = 0.40 GeV has been fixed using the rich experimental data on the 
Bj and i? ± mesons. In the SU(3) limit, this parameter should be the same in B s decays. However, facing the 
high precision experimental data, one has to consider the small SU(3) breaking effect, i.e. the s quark momentum 
fraction should be a little larger than that of the u or d quark in the lighter B mesons. The shape of the distribution 
amplitude is shown in Fig|3]for oj b — 0.45 GeV, 0.5 GeV, and 0.55 GeV. It is easy to see that the larger Wb gives a 
larger momentum fraction to the s quark. We use uji> = 0.50 ± 0.05 GeV in this paper for the B s decays 
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B. B c meson wave function 



For B c meson, we only consider the contribution from the first Lorentz structure, like B q (q — u, d, s) meson, 



/ 2Nc. 



■ m Be )~1h<i>B c (x,b). 



(14) 



Consisting of two heavy quarks (b,c), the B c meson is usually treated as a heavy quarkonium system. In the 
non-relativistic limit, we adopt the model for the distribution amplitude as 29]: 



4>b c (x,b) 



f. 



2J2N, 



S(x — m c /mB c ) exp 



-w 2 b 2 



(15) 



in which exp [— iw 2 6 2 ] represents the kr dependence. /b c and N c = 3 are the decay constant of B c meson and the 
color number respectively. 

For the final state meson wave functions, such as pion and Kaon, we refer the readers to the papers dealing with 
various decay channels 0, [l8j]. For the heavy D meson in the final states, we utilize the heavy quark symmetry to 
simplify the Lorentz structure 



III. NUMERICAL RESULTS AND DISCUSSIONS 

In the numerical calculations, one needs wave function parameters and form factors in the QCD factorization 
approach as input to give numerical results. To deal with the annihilation type diagrams, one needs more parameters 
preferably fitted by experiments. In the soft-collinear effective theory, one needs to fix even more parameters by 
fitting experimental data 0, [l3| , since there are more contributions such as charming penguins without theoretical 
predictions. For the perturbative QCD approach discussed more thoroughly in this paper, one needs only wave 
function parameters. Since some parameters change time to time, one needs to consult the original paper for input 
parameters for each predictions. We do not list them here one by one. 



A. Charmless hadronic two body decays of B s meson 

In the SU(3) symmetry limit, the B s is very similar to the B° meson, which is also called the U-spin symmetry. 
However, the precision of current experimental measurement has already reached the size of SU(3) breaking effect in 
theoretical calculations. In ref.Q, we performed a systematic study of all the charmless B s — > PP, PV, VP and VV 
decays (here P and V stand for the light pseudo-scalar and vector mesons, respectively). After our predictions, some 



of the channels are measured by the experimental data 2l| , which are shown in table [U The theoretical errors for these 



entries correspond to the uncertainties in the input hadronic quantities, from the scalc-dcpcndcncc, and the CKM 
matrix elements, respectively. For comparison, we also cite the theoretical estimates of the branching ratios in two 
kinds of QCD factorization framework: QCDF I Q & QCDF II and in SCET [9]. Among them, the B s -> tt+tt" 
is the first channel of annihilation type B s decays. The measured value is well consistent with our pQCD predictions. 
In the soft-collinear effective theory, the annihilation type diagrams are argued to be small and neglected, thus there is 
no prediction for this pure annihilation type decays. Up to now, there is also no SCET calculations for the two vector 
final state decays. As mentioned in the introduction, the annihilation type contributions arc difficult to predict in 
QCD factorization approach, since there is endpoint singularity in these diagrams. A new QCD factorization analysis 
has been performed after the experimental measurement of B s — > 7r + 7r~, which confirms that a large annihilation 
contribution is needed [T^]. This annihilation B s decay is a further confirmation of other annihilation type B° decays, 
such as B° — > Dj K + , which is also well consistent between pQCD theory and experimental measurements 

In addition to the branching ratios, which have large theoretical and exp erimental uncertainties, there is also a 
first time direct CP asymmetry measurement in the B s — > K~n + decay [191 ] 

A CP = 27 ± 8 ±2%. (16) 
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TABLE I: Some of the CP-averaged branching ratios (xlO -6 ) of B s —¥ PP and VV decays obtained in the pQCD approach 
(This work); the errors for these entries correspond to the uncertainties in the input hadronic quantities, from the scale- 
dependence, and the CKM matrix elements, respectively. We have also listed the updated experimental measurements [21I ]. 
For comparison, we also cite the theoretical estimates of the branching ratios in two kinds of QCD factorization framework: 
QCDF I 0] & QCDF II 0], and in SCET @|. 



Modes 


Class 


QCDF I 


QCDF II 


SCET 


This work 


Exp. 


b"^K+tt- 
B° s -> K+K- 

B s — > 7T+7T 


T 
P 
ann 


1n 9+4.5+3.8+0.7+0.8 
iu - -3.9-3.2-1.2-0.7 
99 7+3.5+12.7+2.0+24.1 
zz -'-3.2- 8.4-2.0- 9.1 

n n9A+ ' 003 + - 025 +°' 000 + ai63 
U.UZ4_ Q03-0 .012-0 .000-0 .021 


r O+0.4+0.4 
J -°_0.8-0.5 
9c 9+12.7+12.5 
— 7.2—9.1 

n 90+0.00+0.10 

u - zu -0.00-0.09 


4.9 ±1.2 ±1.3 ±0.3 
18.2 ± 6.7 ± 1.1 ±0.5 


7 ,,+3.2+0.7+0.5 
'•"-2.3-0.7-0.6 
,„ ,,+4.2+7.5+0.7 
lo.D_ 3 2-4.1-0.2 
n C7+0. 16+0. 09+0. 01 
u '°'-0. 13-0. 10-0. 00 


5.4 ±0.6 
24.5 ± 1.8 
0.73 ± 0.14 


B s -> 4>4> 

B s -> K *°K*° 


P 
P 


91 O + 1-1+30.4 
zl -°-l. 1-17.0 
q +0.5+11.3 
-0.4-6.8 


7+2.6+II.3 
lu - ' -2.1-8.8 
gg+1.1+1.9 
°-°-1.4-1.7 




ntr 0+8. 3+16. 7+0.0 
• J0 - ,J -6. 9-10. 2-0.0 
7 O+1.9+3.8+0.0 
' ' -1. 5-2. 2-0.0 


19 ±5 
28.1 ± 7.2 



Our calculations in ref.[7j give the direct CP asymmetry as 

A CP = 24A±li±litli%, (17) 

while the QCD factorization approach gives a result with a minus sign if not fixing the annihilation digram contri- 
bution |4j 

A _ a 7+2.1+3.1+0.2+15.50/ f-\a\ 
A C P - — D. 2.2 — 2.9 — 0.4— 15.2 '°' U 8 J 

It is easy to see that the pQCD results agree with the experimental measurement quite well, which means that only 
the pQCD approach gives the right sign of strong phase, since the direct CP asymmetry is proportional to the sine 
of strong phase difference. This is a further example of the right prediction of direct CP asymmetry in pQCD after 
the P>° —> K + n~ and P>° —> n + ir~ decays [izoi . The last large theoretical uncertainty in the QCD factorization 
result is from the annihilation type diagram contribution, that is only a non calculable parameter. The later QCD 
factorization with fixed large contribution from annihilation diagrams gives similar results with the pQCD predictions 
|: A C p = 20.7™ 9 8 %. 



B. Hadronic two body B decays with one tensor meson in the final states 



Recently, several experimental measurements about charmless B decay modes involving a light tensor meson (T) 
in the final states have been observed 22|. These decays have been studied in the naive factorization approach 
[23I and relativistic quark model 24 1, with which it can be easily shown that (0 \ \ T) — 0, where j M is the 
(V ± A) or (S ± P) current 25|. The factorizable amplitude with a tensor meson emitted vanishes. So these 



decays are prohibited in the naive factorization approach. The branching rations predicted in the naive factorization 
approach are too small compared with the experimental results, which implies the importance of nonfactorizable 
and annihilation type contributions. The recent QCD factorization (QCDF) approach analysis [3] proved this. It is 
worth of mentioning that the perturbative QCD (PQCD) approach [5( is almost the only method to calculate these 
kinds of diagrams, without fitting the experiments. 

The numerical results of B — > PT decays with AS* = 1, together with Isgur-Scora-Grinstein-Wise II (ISGW2) 
model 23 1 and QCDF results [25 1 are shown in table [IT] The experimental data are from Ref . j/^ . The results 
of B — > PT decays with AS — and also the CP asymmetry parameters for all of these decays can be found 
in ref. 26]. Among the considered B — > PT decays, the PQCD predictions for the CP-averaged branching ratios 
vary in the range of 10 -5 to 10 -8 . From the numerical results, we can see that the predicted branching ratios of 
penguin-dominated B — > PT decays in PQCD are larger than those of naive factorization ;23j by one or two orders of 
magnitude, but are close to the QCDF predictions (25j. For illustration, we classify these decays by their dominant 
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TABLE II: The PQCD predictions of CP-averaged branchin g ra tios (in units of 1CT 6 ) for B -> PT decays with AS = 1, 
together with Isgur-Scora-Grinstein-Wise II (ISGW2) model [23J and QCDF results [25| . The experimental data are from 
Ref.Q. 



Decay Modes class This Work ISGW2 [24] QCDF [4] Expt. 
B+^K?n+ PA 0.9 l^™!^ ~ 



B+ - 


-»■ K* + 7V° 


PA 


n A +0.1 +0.1 +0.1 
-0.0 -0.1 -0.1 


0.090 


— 1.9 




B+ - 




T,PA 


o -i +0.7 +0.6 +0.6 
£ - 1 -0.6 -0.5 -0.5 


0.31 


4.915;* 


< 45 


B+ - 




PA 


q 1 +0.9 +0.9 +1.1 
°- 1 -0.8 -0.8 -0.9 


0.011 


S.41^ 1 




B+ - 


■+ fiK+ 


PA p 11 o +2.7 +3.2 +3.0 
L,CJ\,C 11. _ 2 4 _ 2 .g _2.7 


0.34 


q a +7.8 


1 06+ 028 

l.UU_ 29 


B+ - 


+ f>K+ 


P,PA 


q o +0.4 +0.9 +1.0 
J -° -0.4 -0.8 -0.8 


0.004 


4-011;* 


< 7.7 


B+ - 


-> K*+ V 


PA,P 


n o +0.2 +0.3 +0.3 
V.O _ 2 _o.2 -0.3 


0.031 


6-8lJ 3 7 5 


9.1 ±3.0 


B+ - 


-> K*+ V ' 


PA,P 


19 7 +3.7 +4.5 +4.0 
±z - 1 -3.2 -3.5 -3.5 


1.41 


19 1+20.7 
lz - 1 -12.1 


28.01*; 3 , 


B°- 




PA 


i n +0.2 +0.2 +0.3 
1-U _ .2 -0.2 -0.2 




q q + 8.5 

- J -3.2 


< 6.3 


B°- 




PA 


n R +0.2 +0.1 +0.2 
u ' u -0.1 -0.1 -0.1 


0.084 


1 2+ 4 ' 3 


< 4.0 


B°- 




T,PA 


c n + 1 ' 6 + 1 ' 4 + 1 ' 3 

°- u -1.4 -1.1 -1.0 


0.58 


a 7+ 17 ' 2 

9 - ' -8.1 




B°- 




PA 


r, n +0.5 +0.4 +0.6 
z ' u -0.5 -0.4 -0.5 


0.005 


4.2l«; 3 5 




B°- 




PA,P 


Q o +2.0 +2.5 +2.6 
J - z -1.8 -2.1 -2.2 


0.005 


3.411;? 


n 7+1.3 
z -'-1.2 


B°- 




P,PA 


q 7 +0.3 +0.7 +0.9 
' -0.4 -0.8 -0.9 


0.00007 


3.8l-| 




B°- 




PA,P 


i n +0.2 +0.3 +0.3 
1'U _o.2 -0.2 -0.3 


0.029 


6.6l 8 3 7 5 


9.6 ±2.1 


B°- 




PA,P 


1 1 fi +3.6 +4.2 +3.8 
11-U _2.9 -3.1 -3.1 


1.30 


12 4+ 21 ' 3 


13.7l|i 



topologies indicated through the symbols T (color-allowed tree), C (color-suppressed tree), P (penguin emission) and 
PA (penguin annihilation). Although we include also the W annihilation and W exchange diagram contributions, 
none of these channels has dominant contribution from these two topology. This is different from that of QCD 
factorization approach |25|, where a large annihilation type contribution is introduced by hand to explain the large 
experimental data for the penguin annihilation channels. For the theoretical uncertainties in our calculation, we 
estimated three kinds of them: The first errors are caused by the uncertainties of the decay constants of tensor 
mesons. The second errors are from the decay constant (fs = (0.21 ± 0.02) GeV) of B meson and the shape 
parameter (wg = (0.5 ± 0.05) GeV) in the B meson wave function The third errors are estimated from 

the unknown next-to-leading order QCD corrections and the power corrections, characterized by the choice of the 
Aqcd = (0.25 ± 0.05) GeV and the variations of the factorization scales, respectively. One can find that for most 
channels, the size of three kinds of theoretical uncertainties are comparative. 

There are large theoretical uncertainties in any of the individual decay mode calculations. However, we can reduce 
the uncertainties by ratios of decay channels. For example, simple relations among some decay channels are derived 
in the limit of SU(3) flavor symmetry 

B{B° -> K; tt°) ~ B{B+ -> K; + tt°) ~ i B(B° -> K^ir') 



B(B a -> a-K+) B(B+ -> a+K°) 



B(B+ -> a° 2 K+) B(B° -> a° 2 K ) 
One can find from table [TT] that our results basically agree with the relation given above within the errors. 



(19) 
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For decays involving one tensor meson and one heavy meson, we also give predictions with large branching 
ratios 27|. These B decays include the Cabibbo-Kobayashi-Maskawa- favored B decays through b —> c transition, 



and the CKM-suppressed B decays through b — > u transition. Since there are only tree operator contributions, no CP 
asymmetry appears in the standard model for these decays. Again, the factorizable diagrams with a tensor meson 
emitted vanish in the naive factorization. To deal with the large non-factorizable contribution and annihilation type 
contribution, one has to go beyond the naive factorization, to apply the perturbative QCD approach. 

In charmed B decays, there is one more intermediate energy scale, the heavy D meson mass. As a result, another 
expansion series of mo /tub will appear. The factorization is approved at the leading of mo/rriB expansion 
It is also proved factorization in the soft collinear effective theory for this kind of decays [6J. Among those decays 
predicted in ref. 27], there are some pure annihilation type decays, such as B° — > D~ and B s — > Da-i. There 
are currently no experimental measurements for these decays. For the first time, our perturbative QCD approach 
calculations give sizable predictions of branching ratios at the order of 1CP 5 and 10 -6 , respectively [27j, which may 
be measured soon in the experiments. 



C. Two body hadronic decays of B c meson 

From a theoretical point of view [3| , the non-leptonic decays of B c meson are the most complicated decays due to its 
heavy-heavy nature and the participation of strong interaction, which complicate the extraction of parameters in SM, 
but they also provide great opportunities to study the perturbative and nonperturbative QCD, final state interactions 
and heavy quarkonium property, etc. It is well-known that B c meson is a nonrelativistic heavy quarkonium system. 
Thus the two quarks in the B c meson are both at rest and non-relativistic. Since the charm quark in the final state D 
meson is almost at collinear state, a hard gluon is needed to transfer large momentum to the spectator charm quark. 
In the leading order of m c /ms c ~ 0.2 expansion, the factorization theorem is applicable to the B c system similar to 
the situation of B meson [29j. Utilizing the kx factorization instead of collinear factorization, the pQCD approach 
is free of endpoint singularity. Thus the diagrams including factorizable, nonfactorizable and annihilation type, are 
all calculable. For the charmed decays of B meson, it has been demonstrated to be applicable in the leading order 
of the m d /tub expansion 0, • 

The two-body non-leptonic charmless decays B c — > PP, PV/VP, VV can occur through the weak annihilation 
diagrams only. The pQCD predictions .29] for the branching ratios vary in the range of 10~ 6 to 1CP 8 , basically 
agree with the predictions obtained by using the exact SU(3) flavor symmetry. The B c — > K K + and other decays 
with a decay rate at 10 -6 or larger could be measured at the LHC experiment. For B c — > PV/VP, VV decays, the 
branching ratios of AS" = processes are basically larger than those of AS = 1 ones. Such differences are mainly 
induced by the CKM factors involved: V u d ~ 1 for the former decays while V us ~ 0.22 for the latter ones. 

Analogous to B — >• Krf'^ decays, we find Br(B c —> K + n') ~ 10 x Br(B c — s- K + n). This large difference can be 
understood by the destructive and constructive interference between the r\ q and r\ s contribution to the B c — > K + n 
and B c —> K + if decay. Because only tree operators are involved, the CP-violating asymmetries for these considered 
B c decays are absent naturally. For B c — > VV decays, the longitudinal polarization fractions are around 95% to play 
the dominant role, except for B c — > cj)K* + ( f L ~ 86%). 

For charmless decays involving scalar or axial vector final states, some calculations have already been done, such 



as B c — > AP decays, B c — > SP, SV decays, B c — > AV(VA) decays in the perturbative QCD approach [31| etc. 

We calculate the CP averaged branching ratios and direct CP asymmetries for B c — > Dt s )P decays, together with 
results from the light-front quark model (LFQM) [32| and the relativistic constituent quark model (RCQM) 
shown in table IIII1 some of which with large direct CP asymmetry predictions. Generally, our predictions for the 
branching ratios in the tree-dominant B c decays are in good agreements with that of RCQM model. But we have 
much larger branching ratios in the color-suppressed, annihilation diagram dominant B c decays, due to the included 
non-factorizable diagrams and annihilation type diagrams contributions. 

Other similar decay channels B c —> Df s \V and the double charm decays of B c meson are also calculated in ref. 34 1. 



We find that the transverse polarization contributions in some channels, which mainly come from the non-factorizable 
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TABLE III: CP averaged branching ratios and direct CP asymmetries for B c —¥ D^P decays, together with results from 
RCQM0 and LFQIVlQ. 









channels 


Class 


This work 


RCQM 


LFQM 


This work 


RCQM 


B c -»• D%+ 

B c D+TY° 

B c ->■ D°K+ 
B c ->■ D+K° 
B c D+r, 

B c £>+t/ 
B c ->■ £>+7T° 

B c ->■ D+ K° 
B c -»• 

B c £>+»/ 


T 

C,A 
A,P 
A,P 
C,A 
C,A 
C,P 
A,P 
A,P 
A,P 


9fi 7+3.1+6.0+0.8 
' -3.5-5.6-0.6 
n oo+0. 24+0. 55+0. 06 
u '° z -0. 16-0. 41— 0.01 
A7 + 17.2+2.2+5.4 
a ' -°-9. 1-1. 7-3.6 
AK Q+15- 6+0. 3+7.4 
-±U.»_ 12.3-0.3-4.6 
n qo+0. 15+0. 21+0. 03 
u - yz -0. 15-0. 25-0. 00 
n Q1 +0.12+0.16+0.06 
u - yi -0. 10-0. 20-0. 03 
n a -1 +0.04+0.01+0.02 
-0.04-0. 02-0. 02 

1+0.9+0.3+0.3 
^^-0. 6-0. 3-0. 2 

1 7 q+1-7+0.5+3.3 
1 I •"-1.8-0.6-1.2 

m n + 4 - 9 +°- 4 + 6 - 7 
"^•"-s^-o.s-s.s 


22.9 
2.1 
44.5 
49.3 

1.9 


4.3 
0.067 
0.35 

0.087 
0.048 
0.0067 

0.009 
0.0048 


,11 9+4.5+11.1+0.8 
4:1 - z -4.6-7.8-1.2 
9 n+6. 3+1. 4+15.0 
z -°-3.0-0.8-18.8 
ni g+4.9+7.4+1.8 
01.0_ 2 g_ 3 7 _ x 3 

9 n+0. 4+0. 9+0.0 
z - a -0. 2-0. 5-0.0 
A(\ c+0. 0+18. 4+15. 6 
™'°-2 .9-14 .0-13 .5 
ia n +0. 6+4. 6+15. 9 
i41 ' u -l. 5-5.2-11.9 
/IR 7+1-4+6.3+2.5 
1U- ' -1.4-11.8-2.8 
n+6. 9+5. 3+0.0 
Ja - J -7.2-8.0-0.3 
9 o+0. 0+0. 4+1.1 
^•"-0. 1-0. 7-1. 2 
i 1+0.1+0.2+0.7 
L -- L -0. 0-0.2-0.6 


6.5 
-1.9 
-4.6 
-0.8 

13.3 



emission diagrams or annihilation type diagrams, are large. The predicted branching ratios range from very small 
numbers of 0(1O -8 ) up to the largest branching fraction of 0(1O~ 5 ). The theoretical uncertainty study in the 
pQCD approach shows that our numerical results are reliable, which may be tested in the upcoming experimental 
measurements. 

IV. SUMMARY 

The current running of LHCb and other experiments measure more and more hadronic B s and B c decays, which 
require a precision theoretical study of these decays. We summarize the recent progress in theoretical study of two 
body non-leptonic B s and B c decays. Many of the next-to-leading order or even next-to-next-to-leading order a s 
corrections have been performed in various approaches. Some of them are important in the phenomenological study 
and factorization study of hadronic B decays. We also emphasize the importance of the large power corrections 
in the heavy quark expansion, especially the large annihilation type contributions. These power corrections are 
essential in the study of direct CP asymmetry and polarization fraction study of two vector final states. When 
study the final states containing at least one scalar or tensor meson, these contributions may give the dominant 
contribution, especially in the QCD factorization approach. The reason is that the (axial-) vector current or (pseudo-) 
scalar density in the standard model can not produce a scalar or tensor meson from the vacuum. This makes the 
factorizable contribution to these decays negligible. Encouraged by the recent experimental confirmation of some 
B s decays channels, a lot of study of the charmless and charmed B s and B c decays have been performed. Intensive 
study of decays involving a scalar, axial vector or tensor meson in the final states are also in the market. A future 
experimental measurement of these decays will give a test of the various factorization approaches. 
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